Palladinized graphene oxide-MOF induced coupling of Volmer and Heyrovsky mechanisms, for the amplification of the electrocatalytic efficiency of hydrogen evolution reaction

In this study, a nanocomposite of palladium supported graphene oxide (GO)/metal–organic framework (MOF) was prepared using electroless deposition of Pd on GO followed by impregnation method of Pd@GO and MOF. The prepared materials were characterized with various analytical techniques and their applications as HER electrocatalysts were evaluated using cyclic voltammetry (CV), Tafel plots, and turn over frequencies (TOFs). The HER results showed a radical increment of H2 production in the nanocomposite through the Volmer reaction together with Heyrovsky or Tafel mechanism. This disclosed that the addition of Pd@GO/MOF in the electrolytic system possessed better catalytic characteristics with enhanced current density which may open a new way for hydrogen production and storage via HER.

Preparation of MOF and Pd@GO/MOF composite. MOF was prepared by the following hydrothermal procedure 30 . Concisely, 2.5 mmol (0.525 g) of H 3 BTC dissolved in 10 mL of ethanol and then mixed with 10 mL solution of 4.5 mmol (1.087 g) of Cu(NO 3 ) 2 .3H 2 O and distilled water. The mixture was stirred for 30 min and then transferred to a 23 mL Teflon stainless-steel autoclave and sealed to react for 36 h at 120 °C in the thermostatic drying oven. After, the resulted product was filtered, washed with ethanol three times, and then dried at 50 °C overnight. GO; and Pd@GO and Pd@MOF were prepared using a modified Hummers approach 31 and electroless Pd deposition method 28 , respectively.
Pd@GO/MOF composite was synthesised through the impregnation method 28 of directly mixing Pd@GO with MOF. Briefly, 0.1 g of as-dehydrated MOF (at 150 °C for 1 h) was dispersed in 10 mL DMF. In a separate beaker, 0.1 g of Pd@GO (50 wt.% of Pd@GO loading in the composite) was dispersed in 1.4 mL DMF. The two mixtures were mixed for 24 h at 50 °C. The mixture was then filtered to obtain the desired product which was washed several times with ethanol. The resulted clean product was dried in an oven at 50 °C for 12 h. Materials characterization. Spectrum II spectrometer (PerkinElmer, South Africa) was used to record FTIR spectra from 400 and 4000 cm −1 at room temperature. The structural phases of the materials were determined using XRD (Philips PW 1830, Industrial Analytical (Pty) Ltd (South Africa), CuKα radiation, λ = 1.5406 Å). Thermal stability was evaluated using a thermogravimetric analyser (STA 4000, PerkinElmer, South Africa). The samples were heated from 30 to 500 °C at a rate of 20 °C min −1 . Auriga field emission scanning electron microscope/energy dispersive spectroscopy (FESEM/EDS) from Carl Zeiss (Pty) Limited, South Africa, was used for morphological characterizations and elemental analysis. FEI Tecnai G 2 20 transmission electron microscopy together with energy dispersive x-ray spectroscopy (TEM/EDX) and high-resolution transmission electron microscopy together with selected area electron diffraction (HRTEM/SAED) from FEI Company (America) were used to study the internal morphology at 200 kV accelerating voltage with camera length of 100 mm. X-ray photoelectron spectroscopy (XPS) analysis was conducted in a Thermo Scientific ESCALAB 250Xi (Waltham, Massachusetts, U.S) with a monochromatic Al K α (1487 eV) x-ray beam. The calibration or peak shift analysis was done using the valence band spectra instead of the carbon peak shift since carbon was monitored as part of the analysis. EPSILON

Results and discussion
Structural characterization. The XRD patterns of MOF, GO, Pd@MOF, Pd@GO, and Pd@GO/MOF nanocomposite are presented in Fig. 1. The pattern of MOF ( Fig. 1) is characterized by an intense peak at 2 θ = 12°, which is indexed as the (222) lattice plane with an interplanar d spacing of 7.68 Å (estimated using Bragg's formula: d = (nλ)/(2sinθ), where d = interplanar spacing, n = positive integer, λ = wavelength and θ = scattering angle), typical for HKUST-1 type of MOF and is in agreement with MOF crystal structure information data (CSID) 30 . On the other hand, there was a reduction and shift of peak intensity of MOF phases with the introduction of Pd and the hkl (222) lattice plane with an interplanar d spacing of 7.64 Å. For GO, the pattern is observed to have one broad peak at 2 θ=10.5° (d = 6.70 Å), corresponding to GO (002) crystal face 31,32 . However, it is observed that after the introduction of palladium, there is a shift of this mentioned peak to the right of the pattern, indicating the change in the environment of the oxygen-containing groups as palladium is incorporated on GO (Fig. 1). Characteristic peaks of palladium in Pd@GO are observed at 2 θ of 38 and 45° similar to the one observed on Pd@MOF 28 . It was reported that Pd coated materials have two reflection indexes of (111) and (200) at 2 θ values of 39 and 45°, respectively, as an indicator of Pd cubic phase 28 . The final composite Pd@GO/ MOF consists of phases from both GO and MOF but with decreased and shifted intensities. The hkl indexes (002) and (222), of GO and MOF phases were calculated to have the interplanar d spacing of 8.73 and 7.64 Å, respectively. The palladium peaks are observed with low intensities as an indication of the palladium presence in the nanocomposite. Figure 2 displays the FTIR spectra of MOF, GO, Pd@GO, Pd@MOFand Pd@GO/MOF. The lattice vibrations of as-prepared MOF are in agreement with the literature 33 . In the spectrum of MOF, the coordination of the carboxylate group in the organic linker was observed at ~ 1300 cm −133 . The bands at ~ 1645 and 1450 cm −1 in the MOF spectrum also relate to the asymmetric and symmetric stretching vibrations of the C = O groups in BTC 28,33 . Lastly in the spectrum of MOF, a vibrational mode directly involving the Cu centre and the organic ligand (Cu-O) was observed at ~ 530 cm −120 . The spectrum of GO was reported elsewhere 34 and possessed C-O and C = O vibrations at ~ 1060 and 1735 cm −1 from carboxyl and/or carbonyl groups 34 . Notably, the introduction of palladium onto GO and MOFreveals different interactions. The presence of Pd in Pd@GO is seen by the reduction in band intensities as compared to the GO spectrum 31 . Whereas in Pd@MOF, a new band was observed at ~ 1200 cm −1 , which were corresponding to the characteristic peaks of Pd nanoparticles 4 .
The spectrum of Pd@GO/MOF shows the presence of MOF and Pd@GO functional group as an indication of the Pd supported composite formation.
The thermal gravimetric analyses of the prepared materials are shown in Fig. 3. MOF is thermally stable up to 370 °C. The dehydration of MOF occurs at 100-125 °C and the framework's collapse at 350-370 °C 22 . Graphene oxide was reported to be thermally stable up to 180 °C, this collapse is accredited to the decomposition of oxygen-containing compounds 31,35 . Upon addition Pd on the surface of GO and MOF, the stability was enhanced as presented in Fig. 3 owing to the Pd interacting with GO and MOF, and hence limiting/reducing the loss in those functional groups. The addition of Pd@GO on the MOF surface displayed an improvement of www.nature.com/scientificreports/ stability as related to MOF and Pd@GO. Furthermore, the transitions of both MOF and GO are observed in the final composite, which confirms the incorporation between the two materials. The Differential scanning calorimetry (DSC) results ( Fig. 4) correlates well with the TGA analysis. There is an exothermic peak at 100-125 °C in MOF, which is due to the heat released when the dehydration process occurs 22 . Moreover, an endothermic peak at 350-370 °C is due to the absorption of heat as the frameworks of the MOF material collapses. With the introduction of palladium onto the MOF materials, there is an enhancement in thermal stability of the frameworks 28,31 . The Pd@GO shows an endothermic peak at 210 °C which is similar to the one of GO, due to the absorption of heat as material loses all the functional groups 31,35 . The final composite shows peaks corresponding to the loss of the two-parent materials (GO and MOF). This also confirms the formation of the composite Pd@GO/MOF.    www.nature.com/scientificreports/ in Fig. 5a, MOF demonstrates the characteristic irregular shape of octahedral crystals 22,28 . The inset in Fig. 5a shows that MOF has a smooth surface. On the other hand, the EDS spectrum of MOF (Fig. 5g) discloses the amount of carbon, oxygen, and copper atoms present in the framework, and elemental analyses are given in Table 1. From the SEM image of Pd@GO crystals (Fig. 5c), it is clear to see well-dispersed nanoparticles of Pd on the GO surface without aggregation as compared to the SEM image of GO ( Fig. 5b) 31 This was further confirmed by the EDS spectrum of Pd@GO ( Fig. 5h) also showing the presence of Pd atoms which are dispersed through the GO sheets and Pd atomic percentage is given in Table 1. It can be seen that the Pd nanoparticles were uniformly distributed on the surface of MOF (Fig. 5d). From the image of the Pd@GO/MOF composite, the rough surface due to the introduction of Pd@GO on MOF is observed (Fig. 5e). The high magnification image in Fig. 5f for Pd@GO/MOF composite displays the interconnecting of particles on the surfaces of the MOF structure. This remark is ascribed to the morphology of Pd@GO on the MOF structure as compared to neat MOF. EDS (Fig. 5i) showed the presence of Pd nanoparticles in the composite at 3.23 atomic percentage (Table 1). Figure 6a-e illustrates the TEM images of MOF, GO, Pd@GO, Pd@MOF, and Pd@GO/MOF, respectively, which shows a clear distinction among the texture of the samples. The image of MOF (Fig. 6a) shows the welldefined octahedral shapes of this crystalline material. After the introduction of Pd, the structure of the MOF develops rough surfaces which can be due to the interaction of Pd on the material surface (Fig. 6d) 28 . As seen in  www.nature.com/scientificreports/ Fig. 6c for Pd@GO with respect to the TEM image of GO ( Fig. 6b) 31 , the image shows well-exfoliated graphene nanosheets and wrinkled transparent sheet-like structure. In the structure of Pd@GO/MOF composites (Fig. 6e), the layers have been established as the interchange between Pd@GO sheets and MOF blocks. It was seen that the oxygen groups of GO in Pd@GO link with the metals of copper dimmers in MOF structure to form a composite 26 with reference to the neat MOF ( Fig. 6a) and Pd@MOF (Fig. 6d 28 . It is interesting to observe that again in the composite, the particles conserved the profile of MOF crystals (see high magnification image in Fig. 6f) signifying that the constrain effects caused by layers of Pd@GO resulting in maintaining the shapes of the original MOF 26 . As shown in Fig. 6g-i, the EDX spectra and Table 1 revealed all the elements present in MOF, Pd@GO, and Pd@GO/MOF confirming the incorporation of Pd@GO onto MOF, respectively. It was observed that the presence of small peaks of Na, Mn, and S are due to impurities trapped in the GO and Pd@GO/MOF during the synthesis process as seen in EDX spectra. The presence of unwanted impurities in graphene is known to have a significant impact on their structural properties 31 . Even though some of this novel HER electrocatalysts may be presented as metal-free, the traces of metallic impurities in the composites could play a role in altering their electrochemical properties for hydrogen evolution reaction 23 . In addition, Fig. 7a demonstrates that the phase contrast of the intact crystal of MOF was barely visible. The measured interplanar d spacings of 0.857 nm (8.57 Å, estimated from scale bar measurement) correspond to MOF in (222) orientation in HRTEM. In addition, the crystallinity of the MOF observed in these images shows the clear diffraction spots which are in agreement with XRD discussed above. The HRTEM image (Fig. 7b) of Pd@GO shows an amorphous characteristic of graphene sheets 24 . The Pd@GO/MOF image in Fig. 7c for HRTEM shows the spacing is 1.00 nm (10.0 Å) which is close to the d-spacing of (222) plane of MOF and (002) plane of GO. These results clearly suggest the development of the hybrid structure of Pd@GO/MOF with reference to as-prepared MOF (Fig. 7a,b). www.nature.com/scientificreports/ To understand the elemental composition of Pd@GO/MOF composite and its main constituents (i.e. MOF, GO, and Pd@GO), X-ray photoelectron spectroscopy (XPS) was employed as an analytical tool to identify the valence states of Pd and Cu in the composite materials and also to identify elements that exist within a material to support EDX results. The survey spectrum of Pd@GO/MOF composite in Fig. 8a clearly shows the existence of Cu 2+ together with oxygen, carbon, nitrogen, and palladium as an indication of composite formation. Table 2 provides elemental identity and quantification of the expected elements in MOF, GO, Pd@GO, and Pd@GO/ MOF materials. The presence of Pd with a 0.2 atomic percentage was observed in both Pd@GO and Pd@GO/ MOF composite. The decrease in the amount of Cu from 9.6 at.% in the MOF to 3.5 at.% in the final composite was due to the incorporation of Pd@GO in the MOF structure. As indicated in Fig. 8b, the Cu2p3 satellite peak is observed at 934.8 eV binding energy and it represents the CuO present in the MOF. The C1s binding energy of the composite is observed in Fig. 8c, the three peaks observed on the C1s spectra corresponds to the C-C at ~ 284.   www.nature.com/scientificreports/ C-C peak proves that more of the carbon atoms were put in a hexagonal structure in MOF due to the introduction of Pd@GO. The Pd3d spectrum of Pd@GO/MOF (Fig. 8d) can be due to oxidation of Pd, shown at 338.8 eV.
Electrochemical characterization. Electrochemical behaviours of MOF, GO, Pd@GO, Pd@MOF, and Pd@GO/MOF composite were investigated using Au electrode in TBAP/DMSO electrolyte (0.1 mol L −1 ) as a supporting electrolyte due to its good stability and electrical conductivity of the system. The voltammograms are represented in Fig. 9a. The oxidation and reduction process of the blank gold electrode was observed at around − 0.50 V which was similar to the one reported in the literature 36 . The voltammogram of MOF displays three oxidation processes attributed to the oxidation copper from Cu/Cu 1+ , Cu 1+ /Cu 2+ , and Cu 2+ /Cu 3+ between 0 and 1 V 37,38 . On the other hand, the CV of Pd@MOF in the same condition shows enhancement of both oxidation and reduction peaks as compared to MOF and this clearly indicates the introduction of Pd on the MOF surface. The CV of GO material shows characteristic cathodic reduction peaks at the region of − 0.500 and − 0.750 V (vs Ag/AgCl) due to the electrochemical reduction of epoxy, peroxy, and aldehyde functional groups 31 . The quasi-reversible processes were observed in the voltammogram of Pd@GO in the same potential window due to the electroactive species of GO 39 and Pd metal 39 which is clearly seen in Fig. 9c. The CV of Pd@GO/MOF on the gold electrode still allows the diffusion of the redox mediator (Cu 2+ /Cu + ) through their layers to the electrode surface 28 with the enhancement of oxidation and reduction peaks at around 0.3 V for conversion of Cu/Cu + and Cu 3+ /Cu 2+ , respectively. The scan rate studies of MOF, Pd@GO, and Pd@GO/MOF composite were performed at various scan rates. The CV curves of MOF, Pd@GO, and Pd@GO/MOF composite are shown in parts b-d of Fig. 9, respectively. As shown in Fig. 9b for MOF, and (d) for Pd@GO/MOF, the peak currents increase with an increase in scan rate due to the generation of Cu ions which are used as catalysts to the HER mechanism 37,40 . Figure 9c demonstrates multiscan voltammograms of Pd@GO which shows well-maintained reversible peaks with the increase in scan rate. From the scan rate studies, it was observed that both anodic (I pa ) and cathodic peak (I pc ) currents are directly proportional to the scan rates from 0.1 to 1 V s −1 . It was seen that the oxidation and reduction couples in MOF and Pd@GO/MOF displayed an electrochemical quasi-reversible process because of the change in peak potential (ΔE p ) and I pa /I pc ratios. However, the unity of I pa /I pc ratios with respect to Cu 1+ /Cu 2+ and its reverse couple; and the logarithm-logarithm of the absolute value of the reductive peak current and scan rate with the slope of 0.5 as shown in Fig. 10a and Table 3 indicate diffusion-controlled characters of the oxidation and reduction processes 41,42 . In addition, the Randles-Sevcik equation was used to determine the diffusion coefficient, D, for a quasi-reversible system as given in Eq. (1) 41,42 .
(1) I p = 2.65 × 10 5 n 3/2 ACD 1/2 (v) 1/2 www.nature.com/scientificreports/ where, n, A, D, C, and v are the number of electrons transferred, electrode area in cm 2 , diffusion coefficient in cm 2 s −1 , the bulk molar concentration of the electroactive species in mol cm −3 , and scan rate in V s −1 , respectively. Figure 10b exhibits the increase in current as a function of the square root of the scan rate, ν 1/2 which is consistent with Eq. (1). The D values as presented in Table 1, were found to be 1.39 × 10 −7 , 4.15 × 10 −7 , 1.08 × 10 −4 , 8.8 × 10 −7 and 7.79 × 10 −7 cm 2 s −1 for MOF, GO 31 , Pd@GO, Pd@MOF 28 and Pd@GO/MOF, respectively. The diffusion coefficient is higher than the one of MOF; a slight increase in the D for Pd@GO/MOF composite is attributed to the improved electrical conductivity caused by the presence of GO in the framework. The decrease in D for the MOF material indicates the repulsive nature of MOF 28 . The introduction of Pd@GO on the MOF surface enhanced the diffusion coefficient of the composite compared to the parent MOF. A similar trend was observed in MOF-based polymer composite 21,22 .
Furthermore, it was seen that instead of the material undergoing only a diffusional process, it can also adsorb on the surface of the electrode 43 . This behaviour can be obtained by directly relating the peak current with the surface coverage (Γ) and the scan rate, where F, R, and T are Fardaic constant, gas constant, and temperature as presented in Eq. (2). Figure 10c shows a linear association of current and scan rate where regression line equation and R 2 are presented to determine the surface coverages (Table 3) for MOF, Pd@GO, and Pd@GO/MOF, and their values were found to be 1.32 × 10 −10 , 1.27 × 10 −9 , and 1.09 × 10 −10 mol cm −2 , respectively, confirming the adsorption of the material on the gold electrode.
Hydrogen studies. The electrochemical hydrogen activities of all prepared samples were measured using H 2 SO 4 as a proton source in 0.1 M TBAP/ DMSO electrolyte using a three-electrode system. Figure 11a presents the plot of current vs potential of the blank gold electrode, MOF, Pd@GO, and Pd@GO/MOF as compared to GO 31 and Pd@MOF 28 . As compared to CV plots in Fig. 11a in absence of a proton source, upon addition of H 2 SO 4 , a new cathodic wave on MOF, GO, Pd@MOF, and Pd@GO/MOF appeared at an onset potential of − 0.5 V with an increase in cathodic current. However, this behaviour was not observed in the case of Pd@ GO, which designates that there was no hydrogen evolution. From this observation, it can be deduced that the   22 . In order to substantiate that the current obtained is a source of the cathodic wave, the concentration-dependent studies for MOF, Pd@MOFand Pd@GO/MOF were carried out and the results are presented in Fig. 11b-d, respectively. Moreover, it can be seen that an increase in H 2 SO 4 concentration results in a further increase in current reading at lower onset potential, thus the new catalytic wave is due to the proton source. Furthermore, a small anodic peak in Fig. 11c, d (see enlargement in Fig. 11d) is attributed to the molecular hydrogen adsorption by the material as observed by Monama et al. 28 . The adsorption of produced molecular hydrogen by the material was due to the presence of Pd particles as an indicator of the improved hydrogen spillover mechanism 28 and also the formation of Pd hydride 41,44 . The Tafel plot was constructed from current density-potential data at various concentrations ranging from 0.033 to 0.450 mol L −1 H 2 SO 4 for MOF, and Pd@GO/MOF nanocomposite whereby the current, i, is measured as a function of overpotential, η. Ramohlola et al. 20,22 , reported that the Tafel slope, b, could also serve as an indicator of either Volmer, Heyrovsky, and Tafel in a multi-step proton transfer process while i 0 measures the performance of an electrocatalyst. In this study, the values of b and i 0 were estimated by linear polarization curves and the results are given in Table 4. In addition, another important parameter that can give insights into the reaction mechanism is the cathodic transfer coefficient (1-α) which can be estimated using the Tafel equation (Eq. 3). Figure 12 and Table 4 show Tafel plots and parameters of Au electrode, MOF, and Pd@GO/MOF composite. It was seen that the blank electrode has a Tafel slope of 200.8 mV dec −1 at 0.300 mol L −1 H 2 SO 4 . At the same acid conditions, MOF and Pd@GO/MOF composite showed lower Tafel slope values as compared to the GO 31 and Pd@MOF 28 . The Tafel slopes of the composite were observed to between 169.9 and 158.1 mV dec −1 at 3.00 and 4.50 mol L −1 , respectively. This suggests that the HER on the prepared electrocatalysts proceeds via Volmer mechanisms, with the adsorption of the proton on the electrocatalyst surface (Volmer step) as the ratedetermining step (Scheme 1). It was reported that the Tafel slopes of 105-150 mV dec −1 in HER describes the Volmer mechanism as a rate-determining method 8 . The results are in good consent with the work reported by Kubisztal et al. 47 when studying the HER behaviour of nickel-based molybdenum composite. It was also shown that a Tafel slope of more than 120 mV dec −1 , originates from either the Volmer rate-determining step or the Heyrovsky rate-determining step with high adsorbed hydrogen atom coverage 48 .
Moreover, the charge-transfer coefficient (α) of 0.5 is known to describe the Volmer reaction or the Volmer reaction together with Heyrovsky or Tafel reactions 8 . As given in Table 4, the α values of MOF and Pd@GO/ www.nature.com/scientificreports/ MOF composite are approaching 0.5 confirming the proposed Volmer mechanism as a limiting step 8,47,49 . It was seen that GO/MOF composite exhibited charge transfer coefficients that are close to 0.5, showing that the ratedetermining steps as Volmer coupled with Tafel or Heyrovsky mechanism 31 . Digraskar and co-workers 50    www.nature.com/scientificreports/ hand, the exchange current density (i 0 ) was estimated from the extrapolation of the Tafel curves 20,21 . The exchange current densities obtained for different electrocatalysts are recorded in Table 4 and these values increase with an increase in H 2 SO 4 concentration. The obtained exchange current densities are higher compared to those reported in other studies 8,17,[45][46][47]52 . From this, it can be concluded that the prepared electrocatalyst possesses fast electron transfer, a large surface area, and favourable HER activities 9 .
In addition, the turnover frequency (TOF) which is defined as the number of hydrogen molecules evolved on an active site in a time period, is also employed to evaluate the catalytic efficiency of the materials 31 . TOF quantifies the specific activity of catalytic centres and can be used to compare various materials 53 . Equation 4 presents the method to determine the TOF values of the material 53 : where j denotes current density, M is the mass percentage of prepared materials, F defines the Faraday's constant, whereas m represents the mass per square centimetre of electrocatalysts estimated from BET surface area of MOF (614. 7 m 2 /g 28 ) and GO (423 m 2 /g 54 ). The TOF values for MOF, Pd@MOF, and Pd@GO/MOF composite are presented in Table 4. It was seen that Pd@GO/MOF synthesised materials resulted in the enhancement of the TOF value as the concentration of H 2 SO 4 increases, for example, the value was obtained to be 7.8 mol H 2 s −1 at the highest H 2 SO 4 concentration (0.45 mol L −1 ). The analysis of the TOF data provides information for unraveling the intrinsic activity of electrocatalysts, as shown in Table 4, Pd@GO/MOF exhibits a slightly higher TOF value of 4.7 mol H 2 s −1 than that for GO, MOF, and Pd@MOF under similar conditions of 0.300 mol L −1 H 2 SO 4 . It is worth mentioning that Pd@GO/MOF composite also showed a higher hydrogen generation rate of 7.8 mol H 2 s −1 than MOF at 0.450 mol L −1 H 2 SO 4 and makes it to be the best HER electrocatalyst in acidic medium. Moreover, Fig. 13a, b present the pH dependence of the current density and potential, respectively, for the experimental data obtained in the present study for MOF, Pd@MOF, and Pd@GO/MOF. The slope of the trendline in Fig. 12a (log(i) vs. pH), signifies the HER reaction order which was determined to be ~ 0.8, suggesting a mixed type mechanism of the Volmer coupled with Heyrovsky or Tafel steps 55 as observed in Tafel parameters. Furthermore, Fig. 13b shows a representation of the Pourbaix diagram, and the slope was seen to be less than 0.5. The observed dependence agrees well with the proposed mechanism of the process above.
Electrochemical impedance spectroscopy (EIS) is one of the electrochemical characterisation techniques that can reveal the electrical resistance of several materials. It is an important technique to disclose the conventional electrochemical measurements, interfacial interaction and to study the electrochemical kinetics of the HER process 56,57 . Figure   www.nature.com/scientificreports/ circuit (EEC), which encompasses Voigt RC elements, involving a solution/ electrolyte resistance (Rs indicated as R1), charge-transfer resistance (Rct signified as R2), and constant phase elements (CPE or Q) as shown by the inset of Fig. 14. It was known that the diameter of the semicircle of the Nyquist plot gives information about charge transfer resistance (Rct) of the interface, and its value is obtained by fitting the EEC 55 . As seen in Fig. 14, the semicircles for all materials are smaller than blank, which could be the reason for enhanced electrochemical properties. It was reported in the literature that the radius of high-frequency semicircle divulges hydrogen absorption and low-frequency semicircle represents the HER kinetics 56 . As observed in Fig. 14, the Rct for Pd@ GO/MOF composite is smaller than the MOF value as shown in Table 5 (i.e., from 249 to 37.7 kΩ), indicating the excellent conductivity of the composites due to the presence of Pd@GO. Rct of materials follows the sequence of blank < MOF < GO < Pd@MOF, Pd@GO < Pd@GO/MOF. Accordingly, a small semicircle describes a high HER activity with a small Rct value for the Pd@GO/MOF composite. In addition, EIS results display a reduced charge transfer resistance and improved conductivity because of the presence of graphene oxide support. In contrast, the constant phase element (CPE or C2) was given to evaluate the exposed active surface area of the electrode. As shown in Table 5, MOF exhibits the largest CPE whereas composite has a 0.457 µF CPE value. Consequently, results advocate that the synthesised Pd@GO/MOF composite is a good candidate for hydrogen evolution reaction.

Conclusions
The Pd@GO/MOF nanocomposite as an efficient electrocatalyst was successfully prepared from MOF and Pd@ GO via impregnation procedure. The presence of GO was very important for the support of Pd atoms through the diffusion process to form Pd@GO/MOF and the addition of defects into GO to expose catalytic sites. The influences of Pd@GO on the structural properties of MOF surface, as well as electrochemical activities for hydrogen evolution reaction, were systematically investigated. The observed promotional roles of incorporated Pd@GO on MOF surface are endorsed to possible synergetic effects between Pd and GO with MOF matrix, resulting in an improvement of HER as part of hydrogen production with enhanced TOF values. Moreover, the b and α values displayed that the HER rate-determining step on the studied MOF and Pd@GO/MOF nanocomposite may be the Volmer mechanism or the Volmer mechanism together with one of other two reactions supported by slopes from the plot of log current vs pH and Pourbaix diagram. In this work, the reported strategy of Pd@GO/MOF nanocomposite which exposes the catalytic active sites of the MOF core, while, maintaining the porous carbon skeleton of GO and the presence of Pd nanoparticles, posed a great potential for more applications in hydrogen production and storage technology thorough hydrogen evolution reaction.